Spin-stretching modes in anisotropic magnets: spin-wave excitations in the 

multiferroic Ba2CoGe207 
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We studied spin excitations of the non-centrosymmetric Ba2CoGe207 in high magnetic fields 
up to 33 T. In the electron spin resonance and far infrared absorption spectra we found several 
spin excitations beyond the two conventional magnon modes expected for such a two-sublattice 
antiferromagnet. We show that a multi-boson spin- wave theory describes these unconventional 
modes, including spin-stretching modes, characterized by oscillating magnetic dipole and quadrupole 
moment. The lack of the inversion symmetry allows each mode to become electric dipole active. We 
expect that the spin-stretching modes can be generally observed in inelastic neutron scattering and 
light absorption experiments in a broad class of ordered S' > 1/2 spin systems with strong single-ion 
anisotropy and/or non-centrosymmetric lattice structure. 

PACS numbers; 75.85.-|-t 75.30.Gw 75.10.-b 76.50.+g 



Magnons are collective spin excitations in crystals with 
long-range magnetic order, often investigated by electro- 
magnetic absorption and neutron scattering experiments. 
Both classical and quantum spin- wave theory of S* = 1/2 
systems predict one magnon branch in the spin-excitation 
spectrum for each spin in the magnetic unit cell [l[ . This 
rule about the number of magnon branches is generally 
accepted and experimentally verified for S > 1/2 spin 
systems as long as the conventional spin- wave theory ap- 
plies, requiring that the lengths (i.e., the absolute values 
of the expectation values) of the spins are preserved in the 
excited states and only their orientations change relative 
to the ground-state configuration Q. However, the pic- 
ture of one magnon mode per spin in the magnetic unit 
cell needed to be surpassed in several /-electron com- 
pounds with complicated quadrupolar ordering, such as 
CeBg [3] and UO2 0. 

Recently, additional spin-wave modes have been ob- 
served by far infrared (FIR) spectroscopy Q and inelas- 
tic neutron scattering (INS) Q in Ba2CoGe207, a sim- 
ple two-sublattice easy-plane antiferromagnet (AF) with 
5 = 3/2 spins 0,[1]- This material has attracted much 
interest owing to its multiferroic ground state where del- 
icate magnetic control of the ferroelectric polarization 
0| and chirality [l3| were realized. Moreover, spin 



waves in Ba2CoGe207 exhibits giant directional dichro- 
ism and natural optical activity at THz f req uencies due 
to the large ac magnetoelectric effect 0, [13] ■ A recent 
numerical diagonalization study on finite spin clusters 
found, besides the two conventional AF modes, addi- 
tional spin resonances with peculiar optical properties 
lol 11 1 . Nevertheless, the understanding of the uncon- 
ventional magnon modes and the coupled dynamics of 
spins and electronic polarization on a fundamental level 
remained an open issue. 

In this Letter, we investigate the spin-wave excitations 
in Ba2CoGe207 over a broad photon energy range com- 
bining electron spin resonance (ESR) and high-resolution 
FIR spectroscopy. The largest magnetic field, 33 T, ap- 
plied in this study drastically changes the antiferromag- 
netic spin configuration for any field direction, in contrast 
to former experiments restricted to i3dc < 12 T. The ori- 
entation of _Bdc and the light polarization relative to the 
main crystallographic axes were systematically varied in 
order to map the field dependence and the selection rules 
of the modes. We derive a multi-boson spin- wave theory 
and show that a large single-ion anisotropy plays a key 
role in the emergence of new magnetic excitations in- 
volving the oscillation of spin length, and that the lack 
of inversion symmetry, a necessary condition for the dc 
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FIG. 1. (color online). Magnetic field dependence of the absorption spectra in Ba2CoGe207 below 2 THz for a representative 
set of light polarizations. The spectra are shifted vertically proportional to the magnitude of the field, Bdc- The distance 
between horizontal grid lines corresponds to 20cm~^ in panels (a) and (c), and 30cm~^ in (b) and (d). The direction of _Bdc 
is indicated in each panel and the spectra for different polarizations and propagation directions (k) of light are distinguished 
by the color. Black triangles and green diamonds represent the position of the resonances determined from the FIR and ESR 
spectra, respectively. The grey lines show the field dependence of the modes obtained in our multi-boson spin-wave approach, 
(d) For Bdc -L [001] in some polarization configurations we observed additional modes (open triangles) that are not explained 
by the theory. Two cases (iJ||[110] and iJ||[001]) are shown here. 



and ac magnetoelectric effects, renders these spin- waves 
electric-dipole active. 

Ba2CoGe207 has a non-centrosymmetric tetragonal 
space group, P42im. The magnetic Co^"*" ions are sur- 
rounded by tetrahedra of oxygens compressed along the 
[001] tetragonal axis. Due to the lack of inversion sym- 
metry a coupling between spins and local polarization 



appears [1^ [iSj . This was observed as a magnetic-order 
induced ferroelectricity in this family of materials in- 
cluding BaaCuGezOr gQ, Ca^rSra-^rCoSisOr Q, and 



BaiMnGeaOr [15| 



Here we study spin-wave resonances of Ba2CoGe207 
on high-quality single crystals [s] in the magnetic phase 
at T = 3.5 K. ESR spectroscopy was performed at 75, 
111, 150, and 222 GHz using solid state oscillators, while 
FIR transmission was measured by Fourier transform 
spectroscopy over the region of 0.15-2 THz (0.6-8 me V) 
with a resolution of 15 GHz. 

An overview of the spectra for representative directions 
of the magnetic field is given in Fig. [1] Two sharp peaks 
are present in the zero-field FIR absorption spectra at 
/ ~ 0.5 and 1 THz, in accordance with former studies 
@, • The first is assigned to the usual optical magnon 



branch gapped by magnetic anisotropy of mostly single- 
ion origin. The second is not captured by conventional 
spin-wave theory, and was shown to respond to both the 
magnetic and electric component of light and termed as 
an electromagnon [5J. For B^c along the tetragonal axis 
as in Fig.[TJa) and (c), the ITHz mode shows a V-shape 
splitting with a double-peak structure on the high-energy 
side. The double-peak structure is clearly visible when 
II [010] [see the blue curves in Fig. [TJc)]. The fre- 
quency of the 0.5 THz mode slightly increases in low fields 
[Fig. (Ha) and (c)], however, it turns back after an avoided 
crossing with the lower branch of the 1 THz resonance at 
Sdc ~ 12 T. 

The rotation of B^c from the tetragonal axis to the 
tetragonal plane affects all the modes drastically (Fig.[T|). 
For Bdc II [100], the 1 THz mode is again split into three 
distinct lines. However, they exhibit only a weak soft- 
ening up to a kink at B^c ~ 16 T, from where the res- 
onance frequencies start to increase quickly [Fig. [2b)]. 
The magnon mode at / ^ 0.5 THz becomes silent with 
increasing field for both polarization directions in the 
Faraday geometry (k||i3dc)- 

An additional low-frequency mode appears in the ESR 
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and FIR spectra when B^c is within the tetragonal plane, 
breaking the fourfold rotoinversion symmetry of the lat- 
tice. This mode corresponds to the quasi- Goldstone 
mode of an easy-plane AF when B^c = 0. The frequency 
of this mode is not affected measurably by the orientation 
of i?dc in the plane and follows a linear field dependence 
down to 75 GHz. Hence, the in-plane anisotropy gap is 
less than 75 GHz. 

If i?dc is in the tetragonal plane, the number of ob- 
served resonances exceeds six in some polarization con- 
figurations. Two representative cases are presented in 
Fig.[ljd). The 0.5 THz mode suddenly splits into a sharp 
and a broad feature at B^c = 5 T, while the 1 THz branch 
consists of at least four resonances. At > 12 T, the num- 
ber of modes is reduced. 

As a microscopic model, we consider the Hamiltonian 
below to describe the S = 3/2 spin Co^"*" ions. Follow- 



ing Ref. [llj, we have a large single-ion anisotropy A, but 
we introduce an anisotropic exchange coupling (J and 
Jz) and neglect the Dzyaloshinskii-Moriya term that ap- 
peared to have negligible effect on the excitations: 

T~L — J ^ ^ (^Sf Sj -\- Sf Sj^ -\- Jz ^ ^ Sj -f 

+ gzzhzS^ + g^h^Sf + hySf)] ,(1) 



(id) 



where indicates nearest neighbor pairs, and the x, 
y, and z axes are parallel to the [110], [110], and [001] 
crystallographic directions, respectively, g^x — gyy and 
gzz are the principal values of the g tensor, and ha — 
MB^dc.a are the components of the magnetic field. 

We assume a site-factorized variational wave function 
l^o) = U^eA l^Aii)) n.es ^ describe the long- 

range ordered ground state of two spin sublattices A and 
B [i is the site index). For example, when B^c \\ [HO] 
and Jz/ J ^ 4,, we get a canted Neel state [31 where the 
expectation values of spin components are 



(*xWlSjfl*x(j)) 



3r/(r/ + 1) 



1 



(cos(/?x,sin(^x,0) . (2) 



The two variational parameters rj and ipA = —^b are 
determined from the minimization of the energy {X = 
A,B). We note that f] 1 corresponds to a spin with 
length smaller than 3/2, the consequence of the on-site 
anisotropy. 

This j^I'o) serves as a starting point to study excita- 
tions: we introduce four orthogonal bosons on each site, 
denoted by aj, where = 0, . . . , 3, so that the vari- 
ational ground state is ]\['x(*)) ~ j,j-(i)] vacuum). Any 
product of the operators on a site can be expressed as 
a quadratic form of the four a bosons and they satisfy- 
ing the expected commutation relations. The number of 
bosons on each site is conserved, J2't=o^u x*^!^ x ~ 
and M = 1 for the S = 3/2 spin. The linear flavor- wave 
theory is a 1/M expansion, where the al ^ and al g with 
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FIG. 2. (color online), (a) The energy of the modes for differ- 
ent values of Jz/J in zero field, b- denotes the lu = Gold- 
stone mode. Only the b+ and d+ modes depend on Jz/J- 
The dashed lines indicate modes in the easy-axis AF state 
that forms below A « 2.7J for Jz = 2. (b) and (c) shows the 
imaginary part of the magnetic xTi^i^) ^^'^ electric x||('^) 
dynamic susceptibilities, respectively, for J = Jz- The shad- 
ing above the lines represent the strength of the magnetic and 
electric response. 



v = 1,2,3 play the role of the Holstein-Primakoff bosons 
and describe the excitations, the generalized spin-waves. 



Replacing aj^ and a^j^ with {M - J2t=i"-l,x";y,xV^'^ 
and performing the expansion in 1/M one can follow 
the procedure of the conventional spin-wave theory, and 
we get a Hamiltonian that is quadratic in boson oper- 
ators and straightforward to diagonalize 17|. A sim- 
ilar approach has been used to describe, e.g., CeBg 
[3, the SU(4) Heisenberg modeljlS], the TlCuClg sron 
ladder [3], and multipolar excitations in the spin-l|20j] 
and spin-3/2 [2lj Heisenberg models. 

The spectrum as a function of A/ J in zero magnetic 
field and for zero momentum is shown in Fig. Ufa). It 
consists of six modes, three for each Co spin in the unit 
cell. A finite anisotropy reduces the 0(3) symmetry of 
the Hamiltonian to 0(2), decreasing the number of zero 
energy Goldstone modes from two to one. 

Let us begin from A = 0. Then rj = 1 and the aj 
creates a spin coherent state with maximal spin length 
3/2. In this limit the b± branches correspond to magnons 
of the standard spin- wave theory and they are decoupled 
from the other modes. The c± and d± are local mag- 
netic transitions with ASw = 2 and 3 corresponding 
to Zeemann energies 12J and 18J, respectively, in the 
Weiss field 4 x (3/2) J of the neighboring spins {Sw is the 
spin component parallel to the Weiss field). The c± and 
d± modes are generally silent in neutron, ESR, and FIR 
spectra, as the magnetic dipolar matrix elements vanish 
in the imaginary part of the dynamic magnetic suscep- 
tibility, Imx^™(a;) cx E/ \{f\S^\^o)\^S{Lo - Uf + c^o). 
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FIG. 3. (color online). Motion of the magnetizations (green arrows) and the local electric polarizations (red arrows) in the two 
sublattices (a) for the Goldstone mode (6_ in Fig. [2]) and (b) for the c_ stretching mode. The blue spheres are the oxygens 
forming tetrahedral cages around the central Co ions. The vertical axis is the tetragonal one, while the horizontal axes point 
along [110] and [110]. Apparent tilting of the axes comes from the perspective view. 



These transitions can only be excited by quadrupolar or 
higher order spin operators. 

As we turn on A > 0, 77 increases and the spin length 
decreases in the Neel ground state [Eq. ^]. The modes 
labeled as c_ and c+ in Fig. [2] are spin-stretching modes, 
with spin length oscillating in and out of phase on the 
two sublattices, respectively. Hence c_ is excited by the 
5^ spin operator, with a finite weight in ImxJJ^™(w) that 
vanishes as (A/J)^ when A/J — )■ 0. Most of the weight 
in lmx^^{uj) comes from the low-energy mode, while 
the contribution of to lmx™^{uj) is oc (A/J)'*, so 
that the sum rule /dwlmx™™(w)/a; — g^^/SJ is fulfilled. 
lmx™J^(uj) is zero for all but the Goldstone mode 6_. 

For a large on-site anisotropy (A J,Jz), 77 —> A/3 J 
in the leading order and the = ±3/2 states are sup- 
pressed in the ground state, reducing the spin length to 
1 . We recover the spectra of isolated spins with single- ion 
anisotropy: two modes with energies w/A — >■ and four 
modes with uj — >■ 2A, in agreement with Ref. 111 . 



From the analysis of the dynamic magnetic susceptibil- 
ity it follows that these unconventional spin excitations 
become observable by ESR, FIR, and neutron scatter- 
ing as soon as the single-ion anisotropy gets significant. 
Moreover, if the crystal lattice breaks the inversion sym- 
metry, spin quadrupolar and electric dipole (or electric 
polarization) operators have the same symmetry proper- 
ties. Thus, a new channel opens to excite these modes as 
the electric field of the incident light can directly couple 
to the spin quadrupolar operators, response expressed 
by the Imx^aiui) Ul- Indeed, d_ and c+ modes with 
only magnetic quadrupol moment are excited this way 
and remain silent in Imx™^(a;) irrespective of the A/J 
ratio. The dynamical electric susceptibility hnx'^ai^) 
shows a strong response for most of the modes, as shown 
in Fig. He). 

Our model describes well the magnetic field depen- 
dence of the spin-wave frequencies in Ba2CoGe207, see 
Fig. [1] From the fit of the experimental data we ob- 
tain A = 13.4 K, J = 2.3 K, J^ = 1.8 K, g^^ = 2.1 and 
gxx — gyy = 2.3. Magnetic field larger than 16 T in the 
easy plane is strong enough to drive a transition from 
a canted AF to an almost saturated magnet. This is 
observed as a kink in the 1 THz modes at i?dc ~ 16 T 



[Fig. [UJb)]. The theory also predicts the onset of fully 
saturated phase for _Bdc > 36T applied perpendicular 
to the easy plane (along the tetragonal axis), inducing 
a gap in the Goldstone mode [Fig. HJa)], in agreement 
with Ref. [H. The V-shape splitting of the 1 THz mode 
and the avoided crossing at Bdc ~ 12T for fields par- 
allel to the tetragonal axis is also reproduced correctly 
[see Fig. [IJa) and (c)]. The lowest- lying mode of the 
/ ^ 1 THz branch is theoretically predicted to be weak 
[see dotted grey line in Fig. [2][b)] and does not appear 
in the experimental spectra. The only feature not ex- 
plained by the model is the splitting of the / ^ 0.5 THz 
resonance above _Bdc = 5 T for fields perpendicular to the 
[001] axis. 

The analytical solution of a pure spin Hamiltonian [see 
Eq. dl])] enabled us to fully characterize the excited states 
in terms of spin and polarization dynamics, implying that 
the electric polarization adiabatically follows the sublat- 
tice magnetization vector and does not have its "own" 
dynamics in the energy range of interest. The motion 
of the sublattice magnetization and the local polariza- 
tion in zero field is visualized in Fig. [3] for the Goldstone 
mode and for a spin-stretching mode. The Goldstone 
mode is associated with the oscillation of the polarization 
along the tetragonal axis and has a direct connection with 
the dc magnetoelectric effect. The spin-stretching mode 
shows more complex polarization dynamics [23| . e.g., for 
A — >■ the polarization still oscillates even though the 
magnetic moment is frozen. 

Our theory describes the unconventional spin- wave ex- 
citations in Ba2CoGe207 and provides a guide for spin- 
wave spectroscopy in a broad class of ordered mag- 
nets with strong magnetic anisotropy and/or a non- 
centrosymmetric lattice structure. We expect unconven- 
tional spin excitations to emerge in the dynamic magnetic 
susceptibility whenever a large single-ion anisotropy is 
present in a 5 > 1/2 system. Moreover, if the inversion 
symmetry of the crystal is broken, these new modes can 
have a dielectric response even in the absence of mag- 
netic anisotropy. They should be detected by THz light 
absorption or inelastic neutron scattering via the induced 
magnetic and/or electric dipole moment. 
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SUPPLEMENT 

In this supplement we show in more details the multi- 
boson spin-wave theory. For simplicity, wc consider the 
case when the magnetic field is in the easy plane. For 
magnetic field perpendicular to the easy-plane the ex- 
pressions become rather complicated, and the energy of 
the modes that was shown in Fig. 1 can be obtained nu- 
merically. 



Bosonic representation of the spin operators 

First we introduce the bosons al^ that create the 
= m states of the 5 = 3/2 spin, i.e. |m) = 
I vacuum). The number of bosons on each site is con- 

I 



served, Ylm^ti'^m = ^^'^ M = 1 for the S = 3/2 
spin. Using the four a bosons, the spin operators can all 
be expressed as quadratic forms, for example 



z\2 



3/2 

E 

m=-3/2 
3/2 

m=-3/2 



(3) 
(4) 



S+ = V3 (4/2"l/2 + "-l/2"-3/2) +2^1/20^-1/2(5) 

{S+f = 2VS (ai/2a_i/2 + «V2"-3/2) • (6) 

Next, we apply an SU(4) rotation in the space of al^ 
bosons: 





-a^ - 
— "■0,A — 


: 




V6r?2 + 2 






= (A,A = 


1 




V"! V - 8r? + 2 




— 'A,A — 


1 






v/6??2 + 2 


4 




1 




V - 8r? -h 2 



-3/2 



4/2 + V3ry(e2 



(7a) 



(7c) 



{2v - 1) {ei'^-al,/, - e-i'^-a^,) - VSr, (e^^'^-aL^/^ - e-i'^-a^,)] , (7d) 



and analogous expressions hold for the bosons on the B- 
sites. 



r 



The spins cant in the direction of the field keeping the 
T] parameter unchanged. The limiting cases for r] are 



Variational solution 



In the rotated basis the variational wave functionn cor- 
responds to the I vacuum). The energy per site, as a 
function of r] and ip = ipA = —^b, reads 

E{r],ip) 3(772 + 3) 



N 



4 (3ry2 + 1 

3^(r? + l) 
3772 + 1 



A 



18772(77 + 1)2 
— J cos 2(p 



(37,2 + 1) 
Qxxhx COS(fi. 



(8) 



Minimizing the -^(77, ip) with respect to variational pa- 
rameters ?7 and (p, we get two solutions: (i) canted Neel- 
state, defined via the following set of equations: 



^ _ 3(37; + 1) (772-1) ^ 

3772 + 1 

2477(77 + 1) 
9xxhx = (3^2 _^ X) ^^^'P- 



(9a) 
(9b) 



V= < 



1 + A+o(A2/j2)^ ifA«J; 
[^-1+0{J/A), if A» J. 



(10) 



(ii) For high enough magnetic field, the spins in the A 
and B sublattice become equal and parallel to the field, 
setting if = 0, and 77 is obtained from: 



^ ^ (77 - 1)(377 + 1) 3(377 + 1) (772-1) 

477 ^ 3772 + 1 ^ ' 



77 — )• 1 as the field — 00. 

While the usual procedure is to solve the equations for 
77 and (f, we prefer to express the A as a function of 77 in 
the following. 
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Multiboson spin-waves 



The bx, cx, and dx (i-e. a,^^x, with X — A,B) bosons 
in Eqs. (j7b|) - (|7dp take the role of the Holstein-Primakoff 
bosons. After the 1/M expansion, the spin dipole oper- 
ators are 



^ ^ sm ipx [b'x- Ox 



V3{V - l)i3r, + 1) 
2(3?72 + 1) 



cos fx [cl^ + Cx 



QV ,.377(77+1) . 

^x^M^-^-^smifx 



'M 



S 



X 



iVSy^lr]'^ -477 + 1 
y/iviv + 1) 



^3772 + IVV - 4?7 + 
while the spin quadrupole operators are 



cos ifx (bx - bx 
= {b^X+bx) + 



^/3(77- 1) (377 + 1) 
2(3772 + 1) 

3(77-1)^ 



(4 



2^/3772 + 1a/V - 4?7 + 1 



sni ifix [Cx + Cx 



d^x + dx 



X) 



M 



3 (772 + 3) 2V3r/ 



4 (3r/2 + 1) 37/2 _^ 



677 



■ cos 2if 



X 



?,Tf + 1 

7^/3(7; + 1)(377 - 1) 

V3772 + iVV - 477 + 



= sin 2(px {b\ 



x-bx 



V3(37;2-l) 
3772 + 1 



i2^x ( 



"^x + '^x 



Qiiji — 1)77 



^37/2 + 1^7772 -477 + 1 



sin 2iy9x I — dx 



+ = M, 



677 



3772 + 1 



sin 2ip 



X 



7^3(77 + 1)(377 - 1) 



+ 



67(7; — 1)?7 



= cos2(^x {b\ 



x-bx 



V3 (37/2-1) 

3772 + 1 



sm 



X ^ ^x 



^x^i + ^i^i = VM 



^3772 + 1^7772 -477 

V3(77 + l)(377 -l) 



= cos2(px (c?!: - rfx) 



v/3772 + 1^7772 -477 + 1 



cosiy9x ^&5s: + ^x) + 7\/3sin(px ^c^, 



X '-X 



6(77-1)77 



S-^-S-^ + S'xSl = y/M 



^37/2 + 1^7772 -477 + 1 

V3(77 + l)(377 -l) 



cos (^x ( dx + "^x 



v/3772 + 1^/7772 - 477 
6(7; - 1)77^^ 



= sin ipx + ^x) ~ 7%/3 cos (^x 



X ^x 



V37f + lvV-4^ 



= sin ipx (d^'x + dx 



(12a) 

(12b) 
(12c) 

(12d) 



(12e) 



(12f) 



(12g) 



(12h) 



I 

where terms that are proportional to M and ^/M are The multiboson spin-wave Hamiltonian up to 
shown. quadratic order in bosons reads: 



8 



where Ti^^^ is equal to mean field energy ([8]), T-C^^^ is iden- 
tically zero when ([9|) is satisfied, and the quadratic term 
has the following form for the solution given by Eqs. 



72773(,;+ 1)2 



{■irf + 1) (Iff - 477 + 1) 



J [d\dA + d^ds 



9(77-1)^ 



(37y2 + 1) (7772 _ 4^+1) 

6^377(77+ 1)(77- 1)2 

(3772 + 1) (7772-477 + 1) 



Jz [d\d\ + dAds + d\dB + d^„dA 



Jz {dAbs + dshA + b\dB + b\d\ + fe^d^ + h^eA + d^hg + d^gbAj 



+ 



36\/3772(77 + 1)(77 - 1)2 

(3772 + 1) (7772- 47^+1) '^ l"^"^ + "s"^ + ''^"■^ + "s"^ 



-J (b^d^ + feR^B + 4^-4 + d^nbi 



12772(7/ + 1)^ 



(3772 + 1) [lif -477+1) 



3(7772-477 + 1) 
J. Jcos2<, 



b\h^g + bAbB + b\bB + 



_^ 3(377 + 1)(t/ — 1) — 7^ l!LLljj sin 2.^ (?^acb + ^AC^ + ^^ca + fo^c^ - fo^cs - 6^c^ - 6bca - &bc^ 

(3772 + 1) ^ 



3(377 + 1)2(77 - 1)^ 

(37/2 + 1)2 



J cos2(^ (c^c^ + cacb + cj^cs + c^c^^ + 6J(77 + 1) {c\ca + c^c^^ . 



(14) 



We note that in zero field ((^ — ±7r/2) and parallel spins for the uniform state in high fields, where the variational 
the Hamiltonian separates into two parts, one involving parameters are given by Eq. (HIl) we get 7^(2) = -^(2) _|_ 
b and d bosons, the other only the c bosons. Similarly, -^(2) ^-^^-^ 



-W(2) _ 



and 



(37/2 + 1) (7772 -477+1) 

6V3J,r/(r/ + 1)^^65(7/ -1)2 



(37/2 + 1) (7772- 47/ +1) 



(^dsdA + d\dB + d\d^g + d^d^ 

dAbB + dsbA + b\dB + fo]4ds + b'^B'^A + ^3^^ + d\hB + d^^A 



(37/2 + 1) (7772 -477 + 1) 

l2J,rf{T^ + lf 3 J (77/2 - 4?7 + 1 



(3772 + 1) (7772 -477+1) 3772 + 1 

12J^7/2(77 + l)2 3 J (77/2 - 47/ + 1 



37/2 1 



,fe^(77+l) (9773 -5772 -77+1) _ 6 J(77 + 1)2 (9773 - 57/2 - 77 + 1) \ / . . 

277 (7772- 47/ +1) (3772 + 1) (7772- 47/ +1) )VA-4+B 



3\/3g^xfe.(r/ - 1)^77 36V3 J(7/ - 1)2772(7/ + 1) ' 
(3772 + 1) (7772 -477 + 1) 



7772 -477 + 1 
GgxxhxV^ir] + 1) 



72j7/3(77+ 1)2 



7772-477+1 (3772 + 1) (77/2 - 47/ + 1) 



(bl^dA + ^'V-B + dl^bA + d^6s 



(^d^d^ + d^ds^ 



= + - 6 J(,7 + 1)) (ct + ct c.) 

3J(377 + 1)2(77 - 1)2 / t , , ,t ,t , ,t 



(15) 



(16) 



9 



The spectrum in zero field 

The finite anisotropy reduces the symmetry down to 

0(2), and the Goldstonc mode associated with turning 
the order parameter in the xy plane is desribed by the 



7 



= 7 oc 2?7(?7 + 1) (b\ + h^B + bA + be) 

-^/3(r7-l)2(4 + 4 + rf^ + dB) (17) 

so that [n^^\j^ = [n<-^\j] = o. 



The energy of the remaining 5 modes: In the mode 
the spin moves in the xy plane and is of the same sym- 
metry character as the Goldstone mode, with energy 



^'i- = J- (18) 



The energy of the 6+ and d+ modes that depend on the 
jz = Jz / J is given by the 



J4 



677 + 6 



- j (9??^ + 977^ + 177?2 + 7?7 + 2 - 877^ j,) - 7277^ 



677 + 6 



377^ 



1 



(477' - v^jz + rfjz - 3vjz - jz)m 



equation. The four aforementioned modes consist of b 
and d bosons. 



Finally, the two "stretching" modes that involve the c 
bosons only: 



In the limit of small A the energies are 



+ 1) (9??^ + I8774 - 6773 + 4772 + 577 + 2) 

Wc_ — „ „ ^ J 



3772 + 1 



QVviv + 1) (V + 18/72 + 877-3) ^ 

= vn ^ 



(20) 
(21) 



0, 

24J [3(J - J,) + A] + 0(A2), 
12J + A+^^+0(AVJ'), 

o J 

12J + A-^^+0(AVn 

18J + 0(AVJ^), 
18J + 0(AVJ^). 



(22a) 
(22b) 

(22c) 

(22d) 

(22e) 
(22f) 



The u)bj^ shows the typical square-root behaviour of the 
anisotropy gap on the exchange anysotropy J — Jz and 
single-ion anisotropy A for small gap. 

In the limit of large single-ion anisotropy (A » J,Jz)'- 

tOb_ = 0, (23a) 
Wfc^ = 32J(4J -Jz + 12J2/A) + 0(1/A2), (23b) 

Wc_ = 2A + 7J + 0{J^/A), (23c) 

Wc+ = 2A + J + 0(j2/A), (23d) 

Wd_ = 2A + J + 0( JVA), (23e) 

ujd+ = 2A + 7J + 0( JVA). (23f) 



